EXPONENTIAL DECAY OF DISPERSION MANAGED SOLITONS 
FOR VANISHING AVERAGE DISPERSION 



M. BURAK ERDOGAN, DIRK HUNDERTMARK, AND YOUNG- RAN LEE 

Abstract. We show that any L 2 solution of the Gabitov-Turitsyn equation de- 
scribing dispersion managed solitons decay exponentially in space and frequency 
domains. This confirms in the affirmative Lushnikov's conjecture of exponential 
decay of dispersion managed solitons. 



1. Introduction 

Consider the one-dimensional non-linear Schrodinger equation (NLS) with period- 
ically varying dispersion coefficient 

iu t + d(t)u xx + \u\ 2 u = 0. (1.1) 

It describes the amplitude of a signal transmitted via amplitude modulation of a 
carrier wave through a fiber-optical cable where the dispersion is varied periodically 
along the fiber, see, e.g., [3l [271 EQ] ■ In (jl.ip t corresponds to the distance along the 
fiber, x denotes the (retarded) time, and d(t) is the dispersion along the waveguide 
which, for practical purposes, one can assume to be piecewise constant. 

In fiber optic cables, the information can be transmitted using localized soliton 
pulses in allocated time slots; the presence of a pulse corresponds to "1" and the ab- 
sence of a pulse corresponds to "0" in binary format. Solitary solutions exist due to 
a delicate balance between the dispersion and nonlinearity. In order to minimize the 
interaction between the individual pulses, one needs to keep the pulses sufficiently far 
apart. A draw-back of solitary solutions for this application is that the soliton solu- 
tions with small support have large L 2 norm unless the dispersion constant is small. 
The technique of dispersion management was invented to overcome this difficulty. The 
idea is to use alternating sections of constant but (nearly) opposite dispersion. This 
introduces a rapidly varying dispersion d(t) with small average dispersion, leading 
to well-localized stable soliton-like pulses changing periodically along the fiber. This 
idea has been enormously fruitful (see, e.g., [201 ID El HH [TTl [THl [23l [25]). Record 
breaking transmission rates had been achieved using this technology [23] which is 
now widely used commercially. 
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To study strong dispersion management regime, it is convenient to write 

d(t) = -do(t/e) + d av , 

Here do(t) is the mean zero part which we assume to be piecewise constant, and d av 
the average dispersion over one period , and e is a small parameter. Rescaling t/e to 
t, the envelope equation takes the form 

iu t + d (t)u xx + ed w u xx + e\u\ 2 u = 0. (1.2) 

Since the full equation (jl.2p is very hard to study, Gabitov and Turitsyn suggested 
to separate the free motion given by the solution of iut + do(t)u xx = in (jl.2p . and 
to average over one period, see [8j [9]. In the cas^H when do(t) = 1 on [—1,0] and 
do(t) = —1 on [0, 1] this yields the following equation for the "averaged" solution v 

iv t + ed w v xx + eQ(v, v, v) = 0, where (1.3) 
Q(vi,v 2 ,v 3 ) := J ^{TrViT^TrV^ds, (1.4) 

and T r = e ird *. In some sense, v is the slowly varying part of the amplitude and the 
varying dispersion is interpreted as a fast background oscillation, justifying formally 
the above averaging procedure. This is similar to Kapitza's treatment of the unstable 
pendulum, see [19]. This averaging procedure yielding (jl.3p is well-supported by 
numerical and theoretical studies, see, for example, [HESS SO], and was rigorously 
justified in [31] in the limit of large local dispersion, i.e., as e — > 0. 

One can find stationary solutions by making the ansatz v(t, x) = e l£ujt f(x) in (jl.3p . 
This yields the time independent equation 

-uf = -d w j xx - Q(f, f, f) (1.5) 

describing stationary soliton-like solutions, the so-called dispersion managed solitons. 
Despite the enormous interest in dispersion managed solitons, there are few rigorous 
results available. One reason for this is that it is a nonlinear and nonlocal equa- 
tion. Existence and smoothness of weak solutions of (|1.5p had first been rigorously 
established in [3l] for positive average dispersion d av > 0. In the case d av = 0, the 
existence was obtained in [15], also see [H] for a simplified proof. Smoothness in the 
case d av = was established in [26] . 

Remark 1.1. By a weak solution we mean / G H 1 in the case d av > 0, or / G L 2 in 

the case d av = 0, such that 

- u(g, f) = d av (g', f) - (g, Q(f, /, /)). (1.6) 

for all g G H 1 . Here {g,f) = J R g(x)f(x)dx is the usual scalar product on L 2 (R). 
By a formal calculation, using the unicity of T r in L 2 , we have 

(9,Q(f,f,f)) = Q(9,f,f,f), 



In fact, our method can be extended to more general dispersion profiles. We will address this 
issue together with the case d av > in a later paper. 
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where 

1 

Q{fuf2,h,h) = J J T r h{x)T r f 2 {x)T r h(x)T r U(x)dxds. (1.7) 

R 

The functional Q(fi, fi-, fi) is well-defined for fj £ L 2 (M) due to Strichartz in- 
equality, see [311 do] • 

The decay of the solutions was first addressed by Lusnikov in [21]. He gave con- 
vincing but non-rigorous arguments that any solution / of (11.50 for d av = satisfies 

f{x) ~ \x\ cos(aox 2 + a\x + a2)e~ b ' z ' as x — > oo (1-8) 

for some suitable choice of real constants a,j and b > 0, see also [22J. In particular, he 
predicted that / and / decay exponentially at infinity. For d av = 0, the first rigorous 
x-space decay bounds were established in [ID], where it was shown that both / and / 
decay faster than any polynomial in the case d av = 0. In particular any weak solution 
is a Schwartz function. 

Our main result confirms Lusnikov's exponential decay prediction: 

Theorem 1.2. Assume that d av = 0. Let f G 1? be a weak solution of (|1.5|) . Then 
there exists ^ > such that 

\m\ <e-"M, \f®\<e-M\, 

where f is the Fourier transform of f. 

We have the following immediate corollary. 

Corollary 1.3. Under the conditions of Theorem li.ffi both f and f are analytic in 
a strip containing the real line. 

Remark 1.4. Weak solutions of 11.51 can be found with the help of a variational 
principle. For d av = it is given by 

P\ '■= sup{Q(/, /,/,/)| H/lll = A} (1.9) 

Note Q(f, /,/,/) = /(! f \e it9 -f{x)\ 4 dxdt = \\e iW -f\\% r4 and e**%f is the space- 

[0,1] m 

time Fourier transform of a measure concentrated on the paraboloid {r = k 2 } CM 2 
with square- integr able density /, 

e itd *f(x) = ^JJ e ixk+itT 5(r - k 2 )f(k)drdk. 

Since, by scaling, P\ = Pi A 2 , the variational problem (11.91) yields the best constant 
in the L 4 Fourier extension estimate 

l|e^/||l4 i]X 4<Pi||/|li 2(R) . (1-10) 

for measure with an L 2 density on the paraboloid. Existence of maximizers for the 
variational problem (jl.9p was established in [15], see also [11] . Thus our Theorem 
11.21 and Corollary 11.31 for any weak solutions of (jl.5j) for vanishing average dispersion 
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show, in particular, strong regularity properties for any maximizer of the Fourier 
extension estimate (jl.lOp . 

The inequality (jl.lOp is, of course, closely related to the one-dimensional Strichartz 
inequality 



for which the sharp constant and the maximizers have been classified in [HEUE], an d 
the Fourier extension problem for the sphere for which the existence of maximizers 
and their properties has recently been discussed in [6]. 

In the proof of Theorem 11.21 the central idea is, as in [21 [15] , to obtain suitable 
exponentially weighted a-priori estimates for the weak solution. In [21 [13], the com- 
mutator of the exponential weight with the Schrodinger operator is easily calculated 
since the operator is local. Variations of Agmon's method work also for relativistic 
Schrodinger operators which are nonlocal. However, in these applications one relies 
on the pointwise decay of the corresponding kernel. In our case, a major difficulty 
arises since our operator Q is nonlocal and the kernel of the free Schrodinger evolution 
has no pointwise decay. We overcome this difficulty by using the multi-linear struc- 
ture of Q and the oscillation of the kernel. This is done in Section [21 where we obtain 
exponentially weighted multi-linear estimates for Q. Multi-linear refinements of the 
Strichartz estimate where first established in [5] and later systematically studied in 
|28j . The results of these two papers focus, however, on the Fourier side and do not 
allow exponential weights. More importantly, we require bounds independent of the 
exponential weights, see Theorems 12.21 and 1 2 . 3 1 b elow . Our bounds are refinements of 
the x-space Strichartz estimates which were developed in [10] and used in conjunction 
with well-known Fourier space Strichartz estimates to prove that any weak solution 
is a Schwartz function. We would also like to note that our proof uses only the fact 
that / £ L 2 and as such does not require any of the previous smoothness or decay 
results. 



We start with two alternative representations of Q inspired by the calculations in 





2. A-PRIORI ESTIMATES FOR Q 



Lemma 2.1. 




(2.1) 



l 




e - i (^-^ 2 +^-^)/(4*) /l(r?l)/2(??2)/ 3 (r? 3 )/4(7?4)5(??1 - m + m - m ) dndt 



R4 



Q(/l,/2,/3,/4) 



(2.2) 



i 




e it ^-^+^-^)/ 1 ( 7?1 )/ 2 ( 772 )/3( r?3 )/ 4 (^)5( r?1 _ m + m - m ) d V dt 



R4 
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Proof. Using the formula 



T t f(x) = J e^-y^/^mdy, 



we get 



T t fi(x)T t f 2 (x)T t f 3 (x)T t f 4 (x) 
1 



to _ r , 2+ ^_ ??4 ) /(2t ) e _, (r/i _ ?y2+% _, )4 ) /( 4 4 ) /i(??i)/2(r?2)/3(%)/4(??4)dr? _ 



From which one obtains (|2.ip by performing the x-integration. 
Similarly, one obtains (12, 2D by using the formula 



To obtain exponential decay of dispersion managed solitons we need the following 
'twisted' dispersion management functionals 

Q^(hi, h 2 , h 3 , h 4 ) := Q^Whu e~ F ^h 2 , e~ F ^h 3 , e - F ^h 4 ) 

Q,Ahl, h 2 , h 3 , h 4 ) := Q(e F ^) hl , e - F ^h 2 , e~ F ^h 3 , e'^h,). 

Here X denotes multiplication by x and P = —id x is the one-dimensional momentum 
operator, and 

F^ £ {x):=fi ] X| H,s>0. (2.3) 
1 + e\x\ 

We have the following theorems which are rather surprising at first sight. They are 
the basis for our proof of exponential decay of dispersion management solitons. 

Theorem 2.2. There exists a constant C such that the bounds 

4 

\Q„, £ {h 1 ,h 2 ,h 3 ,h A )\ <C]]_\\hj\\2, (2-4) 

4 

\Q„, £ (h 1 ,h 2 ,h 3 ,h 4 )\ < CjIlNh (2-5) 

hold for all /i, e > 0. 

Theorem 2.3. There exists a constant C such that if for some l,k £ {1,2,3,4} 
r = dist(supp (hi), supp (/ifc)) > 1 then 

C 4 

\Qn,e{hi,h 2 ,h 3 M)\ < -pTTlNIa ( 2 - 6 ) 

^7=i 
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for all n,e>0. Moreover, if t = dist(supp (hi), supp (%)) > 1 then also 

C 4 
v T j=i 

Remark 2.4. The point of Theorems 12.21 and 12.31 is that the constant in the bounds 
is independent of fx, e > 0. We explicitly allow e = 0, which at first seems to be in 
conflict with the fact that we need e F ^ fi h\ G L 2 . However, in this case we can restrict 
ourselves to compactly supported functions hi and then use the a-priori bound and 
the density of these functions in L 2 . 

Let M be a multiplier in the variables 771 , rj 2 , % , 774 and define the oscillatory func- 
tionals 

Kl I (h 1 ,h2,h 3 ,h 4 ) := 
l 

I e - i (.r^+.3 2 -^)/(4t) M( ^ ) 7^/ l2(%) /^/ l4(r/4)(5(r/1 _ m + m _ m) dr]dt 

R4 

K 2 M (h x ,h 2 ,h 3 ,h^) : = 
1 

' e i^ ^^^+^-^)M(r / )7I^M^(r/2)M^/i4(r/4)<5(r/l - m + 773 - 774) dT/dt 



R4 

Note that by Lemma 12 .11 we can rewrite the twisted functionals as 
Q^e{hi,h 2 ,h 3 ,h 4 ) = —K^ e (hx,h2,h3,h4 : ), 

Q»,e{hi, h 2 ,h 3 ,h 4 ) = ^Kl IfiE (hi,h 2 ,h 3 ,h 4 ), 

where 

M^ £ (n) = e^M-F^eM-P^eM-F^Am) _ 
Note that by the triangle inequality the function Mae is bounded by 1 on the set 
Vi ~ V2 + V3 ~ V4: = f° r an y A*) e — 0. Therefore Theorems 12.21 and 12.31 follow 
immediately from the Propositions 12.51 and 12.61 below. 

Proposition 2.5. Let M := ewPrn-r^-Hfy-rii^o M(rfr, m, Vs, Vi) < 00. Then Kfc, 
n = 1,2, is well-defined for all hj G L 2 (W). Moreover, 

4 

^tji^MMM^ <MjJ ll^ilh (2-8) 

where the implicit constant is independent of M and hj, j = 1,2,3,4. 

Proof. By scaling, we can assume M = 1. Let a{rj) := rj\ — n\ + rj 2 — n\. We write 

1 

1 



\Kh\ < 



R4 



-e Mv)/m dt 
t 



\ M (v)\ n l%( r ?i)l <5 ( 7 ?i - % + % - m)dv 
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<mJ\J ~e- ia W^dt\ f[ {hjivjMrjx -V2 + V3- m)dr). (2.9) 

k 4 o i =1 

Now we divide the t-integral into two pieces t < \a(rj)\ where oscillations will be 
important and t > \a(r])\. More precisely, 

/min(lja(r;)|) _ ia ^y^ 4 
/ dt\ H |fy(rfe)|*fai -m + V3- Vi)dv (2.10) 

R 4 ■ J=1 

1 4 

J jU\h J (r] j )\S(r ]l -V2 + V3-m)dv- (2.11) 



R 4 min(l,|a(r;)|) 

Let us introduce the following functionals, which, for later flexibility, we define in a 
little bit more generality than needed at the moment. For any (measurable) subset 
icR 4 let 

r 4 

h(A) := / min(l, K^ 1 ) ]J \hj( Vj )\S( Vl - m + m - m )dr) (2.12) 
A i= l 

1 4 

h{A):=j j Ifllh^Mm-m + m-ri^dridt. (2.13) 
An{|a(j7)|<t} j=1 

By Fubini-Tonelli 

1 4 

<™> = j \ j IT l^faMm - % + % - = / 2 (M 4 ). 

|a(r?)|<t J = 1 

To estimate (|2.10p we employ the following bound, which follows by the change of 
variable r = 1/t and then an integration by parts. 



b 



-e~ ia ^dt 
t 



< « (2.14) 







Using this one sees 
hence 

\K X M \ </!(]R 4 )+/ 2 (IR 4 ) (2.15) 
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We start to estimate the second term. Since \a(r))\ = 2(771 — 772) (772 — 773) on the set 
f]i ~ V2 + f]3 — ??4 = 0, we can estimate the 7/-integral in I 2 (M. 4 ) by 

4 

f 4 

/ II ^ - ^ + m - m)dv (2.I6) 

|77i— t?2|<v^ 5 

4 



\a(r,)\<t i- 1 



7 = 1 

\V2~V3\<Vi 

The first integral on the right hand side of (|2. 16[) can be bounded by 

J |/7i(??i)||/i2(??2)||/i3(%)l|/i4(m -772 + m)\driidmdm 

\rn-m\<Vi 

< HM2IIM2 J \hi(rn)\\h 2 (ri2)\dr]idrj2 

\vi— V2\<Vi 

<\\h 3 \\ 2 \\h42( J Ihiimtfdmdm) 1 2 ( / I^MI 2 ^^) 1 2 

\vl-V2\<Vi \vi-V2\<Vi 
= 2>/t||^l||2||fc2||2||^3||2||/l4||2, 

where we used Cauchy-Schwarz inequality first in the dr/3 integral, then in the dr)\dr) 2 
integral. The second integral can be estimated similarly. Thus 

l 

^ - .. .......... r i 



/ 2 (M 4 ) < ^l/iillzll^Hall^bll^lb J -jft = 8M2IIM2IIMI2IIM2. 



To estimate ii(R ) we split the 77-integral into two disjoint regions 

At :={77GIR 4 : \m ~ m\ < 1 or \m ~ m\ < 1} 

= {77 G R 4 : (t?! - m \ < 1} U {77 € M 4 : |7? 2 - t? 3 | < 1} 
A 2 : = {77 G M 4 : 1 771 — 772 1 > 1 and |??2 - 7? 3 | > 1}. 

Obviously, ii(IR 4 ) = h(Ai) + h{A 2 ). For h(Ax), we bound the minimum in (f2~T2|) 
by 1 and then estimate the remaining integral as in (|2,16p but now for t = 1. This 
shows 

HM) < \\hl\\ 2 \\h 2 \\ 2 \\h 3 \\ 2 \\h4 2 . 
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On the other hand 

h{A 2 ) < / 7—-T-, 0(771 - 772 + 773 - V4)dr] 

J Hv)] 

< f \hi(m)\\h2{m)\\h3(m)\\Mm-m + m)\ d d d 
~ J \vi - mWm - m\ 

l';i-')2l> 1 

<r II/, 11 11/, 11 11/, 11 f /" IM^)! 2 , , y/ 2 

v 7 |»7i - %rh2 - %| y 

l'7l-')2l> 1 

Ito-toI^ 1 

< ll^ilbll^lbll^slbll^lb- 

This finishes the proof for K\,. The proof for is simpler. Using the inequality 



(2.17) 



1 



e iat dt 



^minClJar 1 ), (2.18) 



one realizes ^ Ii(R 4 ) and then proceeds as in the bound of /1 (1R 4 ) . | 

A refinement of this proposition, when at least two of the functions, say, hj and 
hk, have separated supports, is 

Proposition 2.6. Assume that M := sup„ 1 _„ 2+ „ 3 _ r)4=0 M(?7i, 772, 773, 174) < 00 and 
that there exist l,k 6 {1, 2, 3, 4} u>i£/i r = dist(supp (/ty), supp > 1. Then, 

M 

1^(^1^2,^3,^4)1^^1111^112' n = l,2. (2.19) 
Proof. Again we can assume M = 1. For 4 C R 4 let 

1 4 

1(A) := J I y i e -^W/W^| 1^(^)1^ - % + % - %)*/. (2.20) 
A J = 1 

Let J^ fc := {77 G M 4 : \r)i - rj k \ > r}. Then < I(Jf k ). By symmetry in (^20]) . 

it is enough to consider the cases (l,k) = (1,2) and (I, k) = (1,3). First we consider 
the case (I, k) = (1, 2). 

Recalling the definitions (]2. 12|) and (I2.13p . we can, as in the proof of Proposi- 
tion [23J bound I(JT 2 ) by 

I{Jl,2)<h{Jl,2)+h{Jl,2) (2.21) 

In the integral defining 12(^12), we have t > 1 0.(77) | = 2|?7i — ^72 1 1 ^72 — 773 ] > 2T|r7 2 — 773 1 , 
which implies 1 772 — 773 1 < t/r. This yields 

1 4 
4(^1,2) ~ / 7 y IT IfyfoMm -V2 + V3- m)dr]dt 

|»72— »?3 I <*/ T- • J _ 1 
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< 



1 



T 



IMhlMWNhll^ 



4 2, 



where we obtained the last line as in the estimate of f)2. 16j) with \ft replaced by t/r. 
To estimate Ii(Jf 2 ) ^ 

A\ := {\ m - m\ > r} n {\ m - m \ < 1} (2.22) 
A\ := {\m - m\ >r}n {\ m - m \ > 1}. (2.23) 

Then, obviously, 

h{Jl^) = h{Al) + I x {Al). 
Similar to (|2.17p . we bound Ii^A^) by 

Had < ||/ il || 2 ||/ i2 || 2 ||/ i4 || 2 ( / ^r^3^j2 dmdmdvs) 1/2 

l»71-')2l> T 
l')2-'73l> 1 

< ^H/illbll^lbll^lbll^lb 



For estimating Ii(A[) we bound the minimum in (I2T2D by Kr?)]- 1 / 2 to sec 

f 1 4 

Il04[) ~ 7 k - uaM^ - ^i/s n - m + m - v^dn 



< 



1 /" \hi{r]i)h 2 {r}2)hz{m)hi{r}i - m + 



1% - %| 1/2 
|»7a— *73|<1 



drjidr]2drj; 



'3 



^ IIM2IIM2 /• N^%I, , 

I T72 — f?3 1 < 1 

< ||/ll|| 2 ||/l4||2/ /" |/i 3 (%)| 2 ^2^3\ 1 /V /' |/'2('/2)| 2 ^/2^/:M' ~ 



Vt ^ J \v2-m\ 1/2 ' ^ J \m-m\ l/2 

\m-m\<^ I ^72 — ^73 1 < 1 

^^ll^ilbll^lbll^lbll^lb, 

where in the third inequality we used the Cauchy Schwarz bound with respect to drji 
and in the forth inequality with respect to the measure |ry 2 — Tjs\~ l ' 2 drj2drj^. This 
finishes the proof for K\j. 

Again the proof for is simpler. Using (|2.18p and the separation condition 
\Vi ~V2\ > t for all (771,772) in the support of ^1(771)^2(^2) one sees \K 2 M \ < h(Jl >2 ) 
and then proceeds as in the bound of I\(Jl 2 )- 

Now we prove the case (l,k) = (1,3), that is, we assume that the supports of hi 
and /13 are separated by r. In this case we have \K\ d \ < I(J{ 3 ) and \K^\ < I 1 (J 1 " 3 ). 
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T~/2 T II A 

The triangle inequality yields J[ 3 C J 12 U ^2 3 > as subsets of R , hence 

1(^,3) < /(«/i t / 2 2 ) + < -)= n hm* 

V i= i 

and similarly for I 1 (JJ" 3 ). This finishes the proof of the proposition. | 

3. Proof of exponential decay. 
Let / be a weak solution of the dispersion management equation. Let 

\\f\W := ||e F -W/l|2, (3.1) 

with Ffi e defined in (|2.3j) . The main step in our argument is to show that for some 
positive n, ||/||^ )6 is bounded in e > 0. 

Fix r > 1 and define, for an arbitrary function /, 

/< : = /X[-t/3,t/3]j /< := fX[-T,T], f> '■= T,t] c > f~ ■= f< - /<• 

Lemma 3.1. Let f be a weak solution of the dispersion management equation for 
some uj > with \\f\\ = 1. T/ien 

"ll/>IL« < ||/>H!U + e^||/>||^ + ||/>|U, £ e 2 ^(-^ + + e 3 ^(-L + 

where the implicit constant does not depend on fx, e, and r . 

Proof. Since / is a weak solution of the dispersion management equation for some 
u) > 0, we have 

w(<p, f) = Q(<p, f, f, f), for any cp G L? . 
Using this with ip = e* Ffhe f>, we obtain 

"\\f>\\le = Q(e 2F »' e f>JJJ) 

= Q»,e (e F ^f>, e F ^ f, e F ^ f,e F -f). 

Let h := e F ^f. Then 

w||/i>|| 2 = Qn,e(h>,h,h,h). 
Writing h = h > + h < , and using the multilinearity of Que, we obtain 

^||/i>|| 2 = Q IJ ,,s(h > ,h > ,h > ,h > ) + Q /J , )£ (h > ,h < ,h < ,h < ) (3.2) 
+ Q^e(h > ,h > ,h > ,h < ) + Q fl ,e(h > ,h > ,h < ,h < ) 
+ (similar terms with permutations of the last three entrees) 
Note that by Theorem 12.21 we have 

\Q tJ ,,e(h > ,h > ,h > ,h > )\ < ||/i>|| 4 , 
\Q p ,e(h > ,h > ,h > ,h < )\ < ||fc>fllM- 
To estimate the remaining terms, we will further split one of the h< they contain into 

\Q»,e{ h >' h <' h <' h <)\ ^ \QfJ,A h >i h <i h <> h <)\ + \Q/i,e(h > ,h <1 h <1 h r ^)\ 
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< "^>lll|/i<l| 2 ||^<^ll + ll^>ll \\h<: \\ 2 \\hr. 



IT 



using Theorem 12.21 Theorem 12.31 and the fact that the supports of h > and are 
separated by 2r/3. Similarly, 

\Q IJ ,,e{ h >' h >i h <' h <)\ ^ \Qn,e{ h >> h >> h <> h <)\ + |Q^e(^>>^>,^< 3 ^~)| 



< ^ll/l>ll 2 ll/l<llll/l < ll + ll/l>ll 2 ll/l<llll/^ 

IT 



Similar estimates hold for the permurtations. Using these estimates in (|3.2p . we 
obtain 

^IIM 2 ^ ll / i>ll 4 + IIM 3 IIM + ^ll^>ll 2 ll^<llll^«ll + IIM 2 IIMIIIMI (3-3) 



+ 



l^>llll^<ll 2 ll^<ll + ll^>llll^<ll 2 ll^"- 



Dividing both sides by ||/i>|| and using /i<, < e MT /, < e^ T f^, and 
obtain 



(3.4) 
1, we 



\h>\\ < l|/i>ll 3 + IIMV T + IlKlle^f 4= + 



which finishes the proof. 



+ e 



1 



+ 



Proof of Theorem \1.2\ . Step 1. We will first determine r > 1 and we pick fx so that 
e^ r = 2. We can rewrite the bound from Lemma [3. li as (with the notation v = \\h > ||) 



( w - -7= - CIIZ-II) 1 ' " C » 2 ~ C " 3 < + II/-H ) • 



Step 2. Let G(y) = "%v — Cv 2 — CzA Let u m3iX be the maxima of G on 



(3.5) 




[0^o]U[z. 1 ,oo) = G- 1 ([0,GK)]) 
Figure 1. Graph of G{v) and the trapping region G _1 ([0, G(fo)]). 
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Step 3. Let Vq = v ma , x /2, and pick r > 1 so that 

i) c(^= + \\U\\) <min( W /2,G( ! / )) ) 



ii) ||/>|| < vq/2. 
With this choice, we rewrite (|3.5p as 

G(||/>IU,e) < C(^0), (3-6) 

which is valid for any e > 0. This is depicted in figure [TJ 

Step 4. Note that by ii) above and our choice of fi in step 1, we have 

||/>ILi < ||e^||oo||/>|| < < vq. (3.7) 

Finally since ||/>|| M)E depends continuously on e for e > 0, and (|3.7p . the inequality 
(13. 6p shows that ||/>||^ )£ is in the same connected component of G _1 ([0, G(vo)]), that 
i s l|/>IU,e £ [0) ^o] f° r an £ > 0. This implies by monotone convergence theorem that 

||/>IU,0 = sup ||/> H^e < v Q . 

£>0 

This shows that e^'"'/ G L 2 . With the obvious change of notation, a similar argument 
using Theorems E21 E3] for Q^.e shows that e^ f £ L 2 , for some ft, > 0. Finally, the 
pointwise exponential bounds follows from the one-dimensional Sobolev embedding 
theorem, or simply by the following 

e^\f(x)\ 2 = e^ I ^f s \f(s)\ 2 ds 



X 



(■OO 

<2 / e ^l|/( S )||/'( S )|d S <2||e^'l/||||/ / ||<oo. 



X 



Similarly one gets pointwise exponential decay of /. ■ 

Acknowledgements: It is a pleasure to thank Vadim Zharnitsky for instructive 
discussions on the dispersion management technique. 

B. Erdogan and D. Hundertmark are partially supported by NSF grants DMS- 
0600101 and DMS-0803120, respectively and Y.-R. Lee by the National Research 
Foundation of Korea (NRF)-grant 2009-0064945. D. Hundertmark thanks Max- 
Planck Institute for Physics of Complex Systems in Dresden and the Max-Planck 
Institute for Mathematics in the Sciences in Leipzig for their warm hospitality while 
part of this work was done. 

References 

[1] M. J. Ablowitz, G. Biondini, Multiscale pulse dynamics in communication systems with 
strong dispersion management. Opt. Lett. 23 (1998), 1668-1670. 

[2] S. Agmon, Lectures on exponential decay of solutions of second-order elliptic equations: 
bounds on eigen) 'unctions of N-body Schrodinger operators, Mathematical Notes, 29, 
Princeton University Press, Princeton, NJ, 1982. 

[3] G. P. Agrawal, Nonlinear Fiber Optics. Second Edition (Optics and Photonics), Aca- 
demic Press, San Diego, 1995. 



14 M. B. ERDOGAN, D. HUNDERTMARK, AND Y.-R. LEE 

[4] J. Bennett, N. Bez, A. Carbery, and D. Hundertmark, Heat-flow monotonicity of 

Strichartz norms. Anal. PDE 2 (2009), no. 2, 147-158. 
[5] J. Bourgain, Refinements of Strichartz' inequality and applications to 2D-NLS with 

critical nonlinearity. Internat. Math. Res. Notices 1998, no. 5, 253-283. 
[6] M. Christ and S. Shao, On the extremizers of an adjoint Fourier restriction inequality, 

preprint 2010. 

[7] D. Foschi, Maximizers for the Strichartz inequality. J. Eur. Math. Soc. (JEMS) 9 (2007), 
no. 4, 739-774. 

[8] I. Gabitov and S.K. Turitsyn, Averaged pulse dynamics in a cascaded transmission sys- 
tem with passive dispersion compensation. Opt. Lett. 21 (1996), 327-329. 

[9] I. Gabitov and S.K. Turitsyn, Breathing solitons in optical fiber links. JETP Lett. 63 
(1996) 861. 

[10] D. Hundertmark and Y.-R. Lee, Decay estimates and smoothness for solutions of the 
dispersion managed non-linear Schrddinger equation. Commun. Math. Phys. 286 (2009), 
no. 3, 851-873. 

[11] D. Hundertmark and Y.-R. Lee, On non-local variational problems with lack of com- 
pactness related to non-linear optics, Preprint 2010. 

[12] D. Hundertmark and V. Zharnitsky, On sharp Strichartz inequalities for low dimensions. 
International Mathematics Research Notices, vol. 2006, Article ID 34080, 18 pages, 2006. 
doi:10.1155/IMRN/2006/34080 

[13] Hunziker and I. S. M. Sigal, The quantum N-body problem, J. Math. Phys. 41 (2000), 
no. 6, 3448-3510. 

[14] S. Kumar, A. Hasegawa, Quasi-soliton propagation in dispersion-managed optical fibers. 

Opt. Lett. 22 (1997), 372-374. 
[15] M. Kunze, On a variational problem with lack of compactness related to the Strichartz 

inequality. Calc. Var. Partial Differential Equations 19 (2004), no. 3, 307-336. 
[16] M. Kunze, J. Moeser, and V. Zharnitsky, Ground states for the higher order dispersion 

managed NLS equation in the absence of average dispersion, J. Differential Equations 

209 (2005), 77-100. 

[17] C. Kurtzke, Suppression of fiber nonlinearities by appropriate dispersion management. 
IEEE Phot. Tech. Lett. 5 (1993), 1250-1253. 

[18] T. Lakoba, D. J. Kaup, Shape of the stationary pulse in the strong dispersion manage- 
ment regime, Electron. Lett. 34 (1998), 1124-1125. 

[19] L.D. Landau and E.M. Lifshitz, Course of theoretical physics. Vol. 1. Mechanics. Third 
edition. Pergamon Press, Oxford-New York- Toronto, Ont., 1976. 

[20] C. Lin, H. Kogelnik, L. G. Cohen, Optical pulse equalization and low dispersion trans- 
mission in singlemode fibers in the 1.3-1.1 [im spectral region. Opt. Lett. 5 (1980), 
476-478. 

[21] P. M. Lushnikov, Dispersion-managed soliton in a strong dispersion map limit, Opt. 

Lett. 26 (2001), 1535-1537. 
[22] P. M. Lushnikov, Oscillating tails of dispersion managed soliton, JOSA B 21 (2004), 

1913-1918. 

[23] P. V. Mamyshev, N. A. Mamysheva, Pulseoverlapped dispersion-managed data trans- 
mission and intrachannel four-wave mixing, Opt. Lett. 24 (1999), 1454-1456. 

[24] L.F. Mollcnaucr, A. Grant, X. Liu, X. Wei, C. Xie, and I. Kang, Experimental test 
of dense wavelengthdivision multiplexing using novel, periodic- group- delay complemented 
dispersion compensation and dispersionmanaged solitons. Opt. Lett. 28 (2003), 2043- 
2045. 



EXPONENTIAL DECAY OF DISPERSION MANAGED SOLITONS 



15 



[25] L. F. Mollenauer, P. V. Mamyshev, J. Gripp, M. J. Neubelt, N. Mamysheva, L. Griiner- 
Nielsen and T. Veng, Demonstration of massive wavelength-division multiplexing over 
transoceanic distances by use of dispersionmanaged solitons. Opt. Lett. 25 (1999), 
704-706. 

[26] M. Stanislavova, Regularity of ground state solutions of DMNLS equations. J. Diff. Eq. 
210 (2005), 87-105. 

[27] C. Sulcm and P.-L. Sulem, The non-linear Schrddinger equation. Self-focusing and wave 

collapse. Applied Mathematical Sciences, 139. Springer- Verlag, New York, 1999. 
[28] T. Tao, Multilinear weighted convolution of L 2 - functions, and applications to nonlinear 

dispersive equations. Amer. J. Math. 123 (2001), no. 5, 839-908. 
[29] S. K. Turitsyn, N. J. Doran, J. H. B. Nijhof, V. K. Mczcntsev, T. Schafer, and W. 

Forysiak, in Optical Solitons: Theoretical challenges and industrial perspectives, V.E. 

Zakharov and S. Wabnitz, eds. (Springer Verlag, Berlin, 1999), p. 91. 
[30] S. K. Turitsyn, E. G. Shapiro, S. B. Medvedev, M. P. Fedoruk, and V. K. Mezent- 

sev, Physics and mathematics of dispersion-managed optical solitons, Comptes Rendus 

Physique, Academie des sciences/Editions scientifiques et medicales 4 (2003), 145-161. 
[31] V. Zharnitsky, E. Grenier, C.K.R.T. Jones, and S.K. Turitsyn, Stabilzing effects of 

dispersion management, Physica D 152-153 (2001), 794-817. 

Department of Mathematics, Altgeld Hall, University of Illinois at Urbana- 
Champaign, 1409 W. Green Street, Urbana, IL 61801. 
E-mail address: berdogan@math.uiuc.edu 

Department of Mathematics, Altgeld Hall, and Institute for Condensed Matter 
Theory, University of Illinois at Urbana-Champaign, 1409 W. Green Street, Urbana, 
IL 61801. 

E-mail address: dirk@math . uiuc . edu 

Department of Mathematics, Sogang University, Shinsu-dong 1, Mapo-gu, Seoul, 121- 
742, South Korea. 

E-mail address: younglee@sogang.ac.kr 



